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Abstract 

We study hole, electron and exciton transport in a charge transfer system in the presence of 
underdamped vibrational motion. We analyze the signature of these processes in the linear and 
third-, and fifth-order nonlinear electronic spectra. Calculations are performed with a numerically 
exact hierarchical equations of motion method for an underdamped Brownian oscillator spectral 
density. We find that combining electron, hole and exciton transfer can lead to non-trivial spectra 
with more structure than with excitonic coupling alone. Traces taken during the waiting time of 
a two-dimensional spectrum are dominated by vibrational motion and do not reflect the electron, 
hole, and exciton dynamics directly. We find that the fifth-order nonlinear response is particularly 
sensitive to the charge transfer process. While third-order 2D spectroscopy detects the correlation 
between two coherences, fifth-order 2D spectroscopy (2D population spectroscopy) is here designed 
to detect correlations between the excited states during two different time periods. 
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I. INTRODUCTION 


Charge transfer lies at the heart of the function of biological photo-synthetic light har¬ 
vesting complexes, as well as organic solar cells. An optically created electron hole pair 
(exciton) must be separated into free charges to enable the function of the system. The re¬ 
action center or an interface between two materials is the place where the separation process 
takes place, followed by a process that utilized the charges. In addition to these materials, 
charge transfer is also important for the nucleobases in DNA and for model systems. 

Charge transfer in these systems is thought to be mediated by vibrations.- 1 - 2 - Vibrations 
can provide the energy fluctuations needed to bring exciton and charge transfer states in 
resonance, as well as dissipate the excess energy. A model that is often employed to model 
electron transfer in biological and chemical systems uses a single nuclear reaction coordi¬ 
nate, which is in turn coupled to many other degrees of freedom.- 1 ^ The coupling to the 
environment leads to friction, which can influence the dynamics of the electron transfer 
process. 

Electron transfer can be studied experimentally by nonlinear optical spectropscopy,- - - 
and, in particular, two-dimensional optical spectroscopy.- - — Over the past decade, ultrafast 
two-dimensional spectroscopy 1 ^ has been developed as a tool to study energy transfer in 
complex systems in real time.—— Oscillating signatures were found which were interpreted 
as a signature of electronic coherence in light-harvesting systems,—^ as well as conjugated 
polymers .- 7 Recently, it was realized that underdamped vibrations can also explain part or 
all of the observed oscillating signal.— - — The interplay of delocalized exciton states with 
coupling to vibrations greatly influences the energy transport .— - — 7 Similar effects have been 
observed in artificial J-aggregates.— It is natural to ask the question how charge transfer 
processes are reflected in these spectroscopies, and what information can be obtained from 
them. 

Calculated two-dimensional spectra of excitonic systems coupled to vibrations have been 
reported . 29 ' 30 In general, they exhibit cross peaks which correspond to the vibrational side 
peaks in the linear spectrum. It is not clear how this structure in the spectrum can be 
used to learn something about energy or charge transport. In particular, because under¬ 
damped vibrations play an important role in the charge transfer process, on would like to 
use two-dimensional spectroscopy to study this fundamental process. The ultimate goal is to 
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understand how nature uses electronic and vibronic coherence to enable the function of pho¬ 
tosynthetic systems, and how these principles can be used to guide the design of materials 
for light harvesting applications. 

In this paper, we set out to theoretically study two-dimensional infrared spectra of a sys¬ 
tem where excitons and charge transfer coexist.— The charge transfer process is mediated 
by coupling to an underdamped vibration, which, in turn, is coupled to a dissipative bath. 
We calculate, for the first time, the dynamics of charge separation and two-dimensional op¬ 
tical spectra in the presence of vibrations. Although two-dimensional optical spectra of an 
electron transfer system were calculated in previous work , 34 the properties of spectra in the 
presence of both charge and exciton transfer, which is relevant to real systems, is still an 
open problem, which we address here. Our hierarchical equations of motion approach allows 
us to incorporate strong coupling to vibrations as well as a proper treatment of system bath 
coherence, which is important in the modelling of two-dimensional optical spectra. 

The remainder of this paper is organized as follows. In section |TT] we present the model 
used in the calculations. In section m we discuss the calculated linear and two-dimensional 
optical spectra for a system where charge transfer and exciton transport coexist. In section 
m we present our conclusions. 

II. MODEL 

Because the natural and artificial systems that form the topic of this paper have an 
enormously complicated structure, it is impossible to treat all degrees of freedom quantum 
mechanically in a dynamic model. However, in order to understand fundamental concepts 
such as electronic and vibronic coherence, a quantum mechanical model of the functional 
part of the system is required. The usual way out of this problem is to model the system 
of relevant electronic degrees of freedom quantum mechanically, while the vibrational envi¬ 
ronment is treated as a bath. In our case, the system will be the exciton as well as charge 
transfer states. While it is possible to model the environment using the laws of classical 
physics, this is not good enough for our current purpose. The reason is that one of our 
aims is to assess the role of vibronic coherence, which is a quantum effect that exists as 
quantum coherence between the system and the bath.— In order to treat the bath quantum 
mechanically, we employ the hierarchy of equations of motion method. 
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The hierarchy of equations of motion method was initially developed for a system coupled 
to an overdamped vibration. This case is now well known.— - — Overdamped vibrations 
protect long-lived electronic coherence, when the time scale of the vibrational damping is 
treated properly.—^ - — Less is known about the situation with a main system coupled to 
an underdamped vibration , 31 ' 42 - 4 '' 48 which is our focus here. Underdamped vibrations can 
dynamically bring charge transfer and exciton states into resonance, leading to rates of 
irreversible charge transfer that are impossible without vibronic coupling. 

Although it is possible to model a more general vibrational mode by employing the 
hierarchy in the Wigner picture,— - — we here use the simpler model of a harmonic potential. 
The main electronic system of interest, modeled by a Hamiltonian Hg is coupled to a single 
harmonic mode, which in turn is coupled to a bath of infinitely many harmonic modes. These 
modes lead to damping of the primary vibrational mode. We will choose parameters in such 
a way that the damping is in the underdamped regime, so that the vibration vibrates. The 
spectral density for the coupling of the primary vibration to its environment is chosen to be 
of Ohmic form. Through a transformation, the model can be transformed into an electronic 
system coupled to a bath of infinitely many harmonic vibrations, with an altered spectral 
density . 47 

The Hamiltonian of the model is then 


h = h s + J2 


' P 2 a . 1 2 / W 

--+ -m a uj a {x a - 2 

v 2m a 2 


Vf 


( 1 ) 


where a indexes the bath modes, p, m and x are the momentum, mass and coordinate, 
respectively, of the bath oscillator, c is the strength of the coupling of the bath oscillator to 
the system and V is a system operator. Hg is the system Hamiltonian. The last term in 
the Hamiltonian corrects for the bath-induced renormalization. 

Because the model of linear coupling to a harmonic bath corresponds to Gaussian statis¬ 
tics, all information about the system bath coupling is encoded in the spectral density and 
the temperature of the bath. The spectral density of the Brownian oscillator model (under¬ 
damped vibration) is 


J(lo) = 2h\- 




( 2 ) 


(uJq — id 2 ) 2 + 7 2 (U 2 ’ 

which has a characterictic frequency ojq, damping rate 7 and reorganization energy A. 

To construct the equations of motion, one needs the quantum correlation function L(t), 
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which can be calculated in the standard way from the spectral density with the equation 


i r 00 /3Tiuj 

L(t) — L 2 (t) ~ iLAt) — — / d ujJiuj) (coth cos cut — i sin ut). 

TT JO 2 


(3) 


After performing the integration, e.g. by contour integration, on finds the result to be^£ 


hXu ; 


L^i) = ® c —(7/2 -K)t _ e —(7/2+iC )t 


(4) 


and 


L 2 {t) - 



(5) 

( 6 ) 
(7) 


where ( = \Juiq — y 2 /4 and u k = 2irk//3h. 

Because the correlation function is a sum of exponentials, a hierarchy of equations of 
motion for the reduced density matrix can be derived in the usual way. If we write the 
correlation function as 

( 8 ) 

k 

the hierarchy is given by 


p"(t) = - ( <«'/ + y p n (t) 

+ Y\v,P<] 

k 

+ n k (A k Vp n * + A* k p n X v k ) . (9) 

k 

n k ( n k ) refer to an increase (decrease) of the respective index by one. The only differences 
with the original hierarchy for the Drude Lorentz spectral density are the presence of an 
extra dimension and the different coefficients A and 7 . In order to calculate two-dimensional 
spectra, a separate hierarchy is used for the three times t\, t 2 and t 3 . At the moment of 
interaction with the light, all tiers of the hierarchy are multiplied with the dipole operator, 
which leads to the correct preservation of memory over the external interaction, and the 
resulting elements are copied from t± to t 2 or from t 2 to t 3 . During t\ 1 only a coherence 
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FIG. 1: The states in the system Hamiltonian. Horizontal lines denote the HOMO and LUMO 
levels of the two molecules and circles show electrons. The states labeled e\ and e 2 are exciton 
basis states, those labeled ci and C 2 are charge transfer states. 

between the ground state and the excited state is present, during i 2 populations and coher¬ 
ences in the one-particle manifold are included, while during t$ again only coherences are 
considered. 

Now that we established the Hamiltonian for the system coupled to an underdamped 
vibration, it remains to specify the model Hamiltonian for the electronic system. The 
simplest possible model system which is directly relevant for the study of photosynthetic 
complexes, and which contains charge transfer as well as exciton transfer effects, is an 
electronic dimer.—— In order to describe a system with electron, hole, and exciton transfer, 
we take as our system Hamiltonian 

Hs = E ct (|cl)(cl| + |c 2 )(c 2 |) 

+ ^(l e i)( e 2 | + |e 2 )(ei|) 

+ t h (\e 2 ){ci\ + |ci)(e 2 | + |ei)(c 2 | + |c 2 )(ei|) 

+ *e(|ei)(ci| + |ci)(ei| + |e 2 )(c 2 | + |c 2 )(e 2 |), ( 10 ) 

where \e\) and |e 2 ) are the exciton states, |ci) and |c 2 ) are the charge transfer states, E CT 
is the energy of the charge transfer states, J is the excitonic coupling, th is the hole transfer 
and t e the electron transfer. Note that the Hamiltonian has four excited states. In the 
calculations reported in the following section, we set t h = 0.1 or 0 , t e = 0.1 or 0 and 
J = —0.25 or 0. The energy of the charge transfer states is set to E^t = —0.3, while 
the reorganization energy is A = 2.0, the damping 7 = 0.2 and the inverse temperature 
f3 = 1.5. All these values are scaled to the vibrational frequency, which is set to ujq = 1. 
The parameters can easily be rescaled to the values relevant for real systems. Because for a 
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typical photosynthetic system, the largest J-couplings are in the range of 50-150 cm -1 , the 
parameters used here correspond to vibrational frequencies of several hundred wavenumbers. 
Such vibrations are, indeed, ubiquitous in photosynthetic systems. Our model parameters 
are therefore directly relevant to real systems. Although the reorganization energy chosen 
here is rather large, such values are expected, for example, for DNA bases.— The states that 
appear in the system Hamiltonian are shown in Fig. [0 

The final model parameter that needs to be specified is the way the vibrational bath 
interacts with the electronic system. Although our approach can be applied to more general 
system bath interactions, we here choose the coupling is such a way that the vibrational 
bath affects the electron transfer. Therefore, the system part of the system bath coupling is 

v = l(-|ei>< ei | + + | Cl )< Cl | - |e 2 )(c 2 |). (11) 

The renormalization term which contains V 2 is added to the system Hamiltonian, H' s = 
Hs + AH 2 , with the reorganization energy A = / 0 °° du J(cu)/7 tcu. This term corrects the bath 
induces shift in the system parameters. 

For the transition dipole operator, which couples the exciton states to the ground state, 
we have 

V = MlsXeil + \ei){g\) + /Mb)(e 2 1 + \^){g\) (12) 

In particular, we choose the two transition dipoles to be parallel and of equal length, so 
that we can ignore the vector nature of the transition dipoles and set pi = /i 2 = 1. This 
assumption can easily be relaxed if one is interested in the effect of laser polarization on the 
two-dimensional spectra. 

III. RESULTS AND DISCUSSION 

A. Linear spectra 

First, we turn our attention to the linear spectra. Because population dynamics only 
affects the line shape of the spectrum, this technique is not very sensitive to the details 
of the coupling of exciton to charge transfer states, i.e. of the charge transfer mechanism. 
However, linear spectra already contain a wealth of information. The peak positions reveal 
the energies of underlying states, while their shape reflects the interaction with the environ¬ 
ment. Although charge transfer states, which are optically dark, are not directly visible by 
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FIG. 2: Linear spectrum with (i) hole transfer only, (ii) electron transfer only, (iii) exciton transfer 
only and (iv) all three modes of transfer together, as indicated above the respective panels. 

themselves, they lead to observable effects on the linear spectra because they couple to the 
bright exciton states. In addition to exciton states, underdamped vibrations have a directly 
observable effect on the linear spectrum. They lead to vibrational side bands in addition to 
the main exciton peak, with positions determined by the vibrational frequency and intensi¬ 
ties derived from the Huang Rhys factors. What is not known, and is an important goal of 
our simulations, is the information content of linear spectra in the case where both charge 
transfer and coupling to underdamped vibrations is present. 

In the linear spectra, shown in figure [2] we observe an interesting effect. In order to 
identify the effect of charge transfer, we compare spectra where (i) only hole transport is 
present, where (ii) only electron transport is present, where (iii) only exciton transport is 
present, and where (iv) all mechanisms contribute. At a first glance, we observe that a 
main peak and vibrational side bands are present in all spectra. The spectra for (i) hole 
(th = 0.1, t e = J = 0) or (ii) electron coupling (t e = 0.1, th = J = 0) only look very 
similar, apart for a weak shoulder at the red side of the main peak in the hole transfer 
spectrum. When we consider the other two spectra, however, some differences are observed. 
The spectrum for (iii) exciton transfer (J = —0.25,th = t e — 0) only is somewhat larger in 
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FIG. 3: Population dynamics of exciton states (black and blue line) and charge transfer states (red 
and green line) for (i) hole transfer only, (ii) electron transfer only, (iii) exciton transfer only and 
(iv) the complete Hamiltonian. The initial condition is all population on one of the exciton states. 
Time is in units of l/wo- 

intensity and shifted. Also in this case, the vibrational side bands are similar to the hole and 
electron transfer cases. However, a surprise is present in the case where (iv) all three modes 
of transport are present (th — t e — 0.1, J = —0.25): the shoulder on the red side of the main 
peak is now very pronounced. This shows that electron, hole and exciton transfer cannot be 
considered separately, but must all be included in order to reproduce the full result already 
for the linear spectrum. The observation of an extra peak in the linear spectrum of a charge 
transfer system coupled to an underdamped vibration is the first main finding of this paper. 

B. Dynamics 

We first look at the dynamics for (i) hole transport only, (ii) electron transport only, 
(iii) exciton transport only and (iv) the complete picture with all tranfer modes combined. 
The dynamics can be understood as taking place in a quantum network— coupled to an 
underdamped vibration. In this picture, a clear distinction can be seen between coherent 
hole transport on the one hand and mostly incoherent electron and exciton transport on 


9 






(i) hole 


(ii) electron 



FIG. 4: Two-dimensional spectrum with (i) hole transfer only, (ii) electron transfer only, (iii) 
exciton transfer only and (iv) all three modes of transfer together, as indicated above the respective 
panels. The waiting time is 10.0 / ojq . 

the other hand. Figure [3] shows the dynamics, which follows intuitive behavior. Vibrational 
oscillations are weak, although they are visible in the exciton transfer. Electron and hole 
transport couple the initially populated exciton state to a charge transfer state. Exciton 
transport depopulates the initially excited state and transfers the population to the other 
exciton state. The different dynamics can be explained as follows. Hole transport (i) is 
not affected by the coupling to the vibration, and therefore exhibits coherent Rabi oscilla¬ 
tions. Electron transfer (ii), on the other hand, is damped due to the interaction with the 
vibration, and is incoherent. Exciton transport (iii) is also mostly incoherent, as a result of 
the effect of the damped vibration. When all three transfer mechanisms are present (iv), a 
combination of partial coherent oscillations and partial incoherent transport are observed. 
Because the system is in a population or excited state coherence during the waiting time of 
a two-dimensional experiment, one would expect to see similar behavior when plotting the 
waiting time dependence of peaks in the two-dimensional spectrum. Because the dynamics 
in the four cases investigated here are clearly different, one may hope that two-dimensional 
spectroscopy can distinguish them. We will see, however, that this is not the case. 
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FIG. 5: Main diagonal peak intensity as a function of waiting time for (blue circles) hole transport 
only, (red crosses) electron transport only, (black stars) exciton transport only, and (green line) 
the full Hamiltonian. Time is in units of 1 / ojq . 

C. Two-dimensional spectra 

In figure [4] we present calculated two-dimensional correlation spectra. Because the effect 
of excited state absorption complicates the discussion, and leads to peaks of opposite sign 
which can usually be separated from bleaching and stimulated emission peaks, we focus 
here on the latter two contributions.— The inclusion of excited state absorption is left as a 
possible extension in future work. 

The first obvious effect in the two-dimensional spectra is the presence of the extra feature 
due to the interplay of charge transfer with vibrations observed already in the linear spec¬ 
trum. Here, as shown in figure [4] (d) we observe the same peak appear as a cross peak with 
the main exciton absorption feature, which partially overlaps with this main diagonal peak. 
In addition, we observe that the vibrational cross peak in the spectrum with all interactions 
present is weaker than when only electron, hole, or exciton transfer is present. Similar to 
our analysis of the linear spectrum, we conclude that cross peaks in two-dimensional spectra 
can not be interpret as arising from charge transfer, exciton transfer or vibrations alone, but 
that they are the result of a complex interplay of all these ingredients. 
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FIG. 6: Vibrational cross peak intensity as a function of waiting time for (blue circles) hole trans¬ 
port only, (red crosses) electron transport only, (black stars) exciton transport only, and (green 
line) the full Hamiltonian. Time is in units of \/ uiq . 

D. Waiting time dependence 

In order to analyze the two-dimensional spectrum further, we plot the intensity of the 
diagonal peak ((cui,^) = (0,0)) and one of the vibrational cross peaks ((cui,^) = (1,0)) 
as a function of the waiting time, shown in figures [5] and [6] Because populations of exciton 
and charge transfer states (or, indeed, coherent superpositions of these) are present during 
the waiting time, one would expect the dynamics of these states to be reflected directly in 
these time traces. We would expect to see clear coherent oscillations in the case of hole 
transport, corresponding to the oscillations found in the dynamics, while such oscillations 
should be absent in the mostly incoherent electron and exciton transfer. However, we hardly 
see this effect in the calculated two-dimensional spectra. In fact, all time traces show similar 
oscillations. We conclude that, in contrast to the dynamics, oscillations in the spectra are 
dominated by vibrational coherence. This is the second main hireling of this paper. The 
fact that the electron and hole dynamics are invisible can be attributed to the fact that the 
charge transfer states don’t couple directly to the light, but interact only via the exciton 
states. 
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FIG. 7: Fifth-order response function as a function of the population times t% and for (i) hole 
transport only, (ii) electron transport only, (iii) exciton transport only and (iv) the full Hamiltonian. 
Note the opposite sign of the response in the case of exciton dynamics. Contours in red and blue 
represent positive and negative values, respectively. Coherence times t\ and £5 are set to zero, while 
£3 = 20. The double-sided Feynman diagram corresponding to the response function calculated 
here is shown to the right. 

E. Fifth order response 

We continued our search for experimentally observable signatures of the charge transfer 
process by considering the fifth-order nonlinear response. 5 ' 58 This order is known to be more 
sensitive to anharmonicity in certain cases.— The fifth-order response functions contain three 
coherence times and two population times. Because correlated population dynamics during 
two time intervals is present, one might hope to see more details of the coupling to charge 
transfer states. 

In figure Owe plot a fifth-order response function in the k s = —k\ + k 2 — k 3 + k A + k 5 
phase matched direction. This spectroscopy can be called ”2D population spectroscopy”, 
because it exhibits the correlations between population dynamics during two separate time 
intervals. The difference between the different cases is clearly visible, in fact, the response 
with only exciton transport present has a different sign. Although we calculated only one 


13 




























contribution to the complete fifth-order signal here this hireling indicates that the fifth order 
response is very sensitive to the charge transfer process. This is the third main result of this 
paper. 

Key to this analysis is the fact that changes in population, as observed in figure [3] are 
reflected in the 2D population spectrum. Because of correlations building up during the 
coherence time £ 3 , the signal is not symmetric along the £2 = £4 hire. We note that it is 
important to set £3 not equal to zero for the observation of the charge transfer dynamics. 
While not studied in detail here, it was found in previous work that coherent dynamics 
during the £5 period is also important for the signal.— By changing £ 5 , it is possible to tune 
the interference of signals arising from excited state or ground state dynamics to enhance or 
surpress the oscillatory motions of specific physical processes. Although we fixed £1 = £5 = 0 
here, we can utilize a representation similar to the 2D frequency domain spectrum shown in 
figure [4] to further visualize the fifth-order response in order to analyze the coherence and 
population dynamics of each cross peak separately. We leave a further investigation of the 
fifth-order signal as an interesting direction for future work. 

IV. CONCLUSIONS 

In summary, we have studied a system in which exciton and charge transfer exist together 
with coupling to an underdamped vibration. This work shows that incorporating electron 
and hole transfer into an excitonic model can lead to additional peaks in optical spectra. This 
peak, which is present in the linear spectrum as well as as a cross peak in the two-dimensional 
spectrum, can not be explained by exciton or charge transfer alone, but is the result of a 
complex interplay of both these ingredients with vibrational coupling. This finding shows 
that one has to be careful in the assignment of peaks in experimental spectra to excitonic or 
vibrational peaks alone, while, in fact, charge transfer may also be an important ingredient. 

Oscillations as a function of the waiting time in the two-dimensional spectrum do not 
reflect the population dynamics and are attributed to vibrations, even though the effect of the 
vibration on the dynamics is weak. Although, therefore, extracting information about the 
charge transfer process from two-dimensional optical spectroscopy is difficult, we fold that 
the fifth-order reponse, in the form of 2D population spectroscopy, is particularly sensitive 
to this process. 
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Our results should help the interpretation of spectra of photosynthetic light harvesting 
systems,—^ conjugated polymers,— model dimers^— and DNA.— Further extension of the 
model presented here to calculate the specific properties of these systems in more detail is 
relatively straightforward and is left as a direction for future work. We also hope that our 
results will stimulate experimental and theoretical work on fifth-order nonlinear processes 
as a probe of charge transfer. An investigation of the full fifth order response for a model 
system is under way. 
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